Observations of the December 13 and 14, 2006, Solar Particle Events in the 80 
MeV/n - 3 GeV/n range from space with PAMELA detector. 

O. Adriani/'2 G. C. Barbarino,^''' G. A. Bazilevskaya,^ R. Bellotti,'^''^ M. Boezio,^ 
E.A. Bogomolov,^ L. Bonechi/'^ M. Bongi,^ V. Bonvicini,^ S. Borisov/^'^'^'^^ S.Bottai,^ 
A. Bruno,^'"^ F. Cafagna/ D. Campana,^ R. Carbone,"^'^^ P. Carlson/^ M.Casolino/"'^^ 
. G. Castellini/^ L. Consiglio,^ M. P. De Fascale,^^'^^ C. De Santis/^'^^ N. De Simone,^^'^^ 

^ : V. Di Felice/°'" A. M. Galper/^ L. Grishantseva/^ W. Gillard/^ G. Jerse,^'^^ A. V. Karelin.^^ 

S. V. Koldashov/^S. Y. Krutkov,9A. N. Kvashnm,^ A. Leonov/^, V. Malakhov, ^^ L. Marcelli/° 
A.G. Mayorov/2 W.Menn,!^ V. V. Mikhailov/^ E. Mocchiutti,^ A. Monaco^ N. Mori/'^ 
(SI ; N. Nikonov,9'i°'ii G. Osteria,^ F. Palma/°'ii P. Papini,^ M. Pearce/^ P. Picozza/^'^^ 

^ ■ C. Pizzolotto,* M. Ricci/^ S. B. Ricciarim,^ R. Sarkar,*^ L. Rossetto/^ M. Simon/^ 

R-Sparvoli/*^'^^ P. Spillantini/'-^ Y. I. Stozhkov,^ A. Vacchi,^ E. Vannuccini,^ G. Vasilyev,^ 
P^ ! S. A. Voronov/2 j. Wu/^ Y. T. Yurkin/^ G. Zampa,^ N. Zampa,^ and V. G. Zverev^^ 

^: Abstract 

o 

^. 

C/5 , University of Florence, Department of Physics, 1-50019 Sesto Fiorentino, Florence, Italy 

'— ^ ^INFN, Sezione di Florence, 1-50019 Sesto Fiorentino, Florence, Italy 

^University of Naples "Federico 11", Department of Physics, 1-80126 Naples, Italy 
Q>^ ■ *INFN, Sezione di Naples, 1-80126 Naples, Italy 



in 

o 



% 



^Lebedev Physical Institute, Leninsky Prospekt 53, RU-119991 Moscow, Russia 

^University of Bari, Department of Physics, 1-70126 Bari, Italy 

■^INFN, Sezione di Bari, 1-70126 Bari, Italy 

*INFN, Sezione di Trieste, 1-34149 Trieste, Italy 

^loffe Physical Technical Institute, RU-194021 St. Petersburg, Russia 

^''INFN, Sezione di Rome Tor Vergata, 1-00133 Rome, Italy 

^^ University of Rome Tor Vergata, Department of Physics, 1-00133 Rome, Italy 

^^National Research Nuclear University "MEPhI", RU-115409 Moscow, Russia 

^■^KTH, Department of Physics, and the Oskar Klein Centre for Cosmoparticle Physics, AlbaNova University 
Centre, SE-10691 Stockholm, Sweden 

"IFAC, 1-50019 Sesto Fiorentino, Florence, Italy 

^^University of Trieste, Department of Physics, 1-34147 Trieste, Italy 

^''Universitat Siegen, Department of Physics, D-57068 Siegen, Germany 

^■"INPN, Laboratori Nazionali di Frascati, Via Enrico Fermi 40, 1-00044 Frascati, Italy 

^^RIKEN, Advanced Science Institute, Wako-shi, Japan 



We present the space spectrometer PAMELA observations of proton and helium 
fluxes during the December 13 and 14, 2006 solar particle events. This is the first 
direct measurement of the solar energetic particles in space with a single instrument 
in the energy range from ~ 80 MeV/n up to ~ 3 GeV/n. In the event of December 
13 measured energy spectra of solar protons and helium were compared with results 
obtained by neutron monitors and other detectors. Our measurements show a spectral 
behaviour different from those derived from the neutron monitor network. No satis- 
factory analytical fitting was found for the energy spectra. During the first hours of 
the December 13 event solar energetic particles spectra were close to the exponential 
form demonstrating rather significant temporal evolution. Solar He with energy up to 
1 GeV/n was recorded on December 13. In the event of December 14 energy of solar 
protons reached 600 MeV whereas maximum energy of He was below 100 MeV/n. The 
spectra were slightly bended in the lower energy range and preserved their form during 
the second event. Difference in the particle flux appearance and temporal evolution 
in these two events may argue for a special conditions leading to acceleration of solar 
particles up to relativistic energies. 

Subject headings: space spectrometer, solar particle emission, energy spectrum 



1. Introduction 



The PAMELA spectrometer (jPicozza et al.l 120071 ) is a space instrument designed for the study 
of primary charged particles and antiparticles in a wide energy interval, mainly from tens of MeV to 
~ 1.2 TeV for protons. It was launched in an elliptical orbit at an altitude between 350 and 610 km 
and inclination of 70 degree in June 2006. The main scientiflc go al of PAMELA is the measurement 
of the particle and antiparticle component in galactic cosmic rays (| Adriani et al.ll2009allbl) . the study 



of galactic cosmic ray modul ation by solar activity, and solar energetic particles (jCasolino et al 



2006al : lDe Simone et al.ll2009l ). This paper reports the PAMELA measurements of the solar proton 
and helium fluxes in the energy range from below 100 MeV/n to several GeV/n during the December 
13 and 14, 2006 solar particle events. 

The problem concerning mechanism and site of solar energetic particles (SEP) acceleration 
remains an open question. Certainly, SEP may be produced after powerful explosive events on the 
Sun, accompanied b y solar flares, Coronal Mass Ejections (CME), bursts of solar X /gamma-rays 
and radio emission (JReamesI 11993 ) . It is clear that not a single mechanism is involved into the 
SEP generation. Stochastic acceleration, shock acceleration, and acceleration by the DC electric 
fields in the process of magnetic reconnection are the main candidates. Acceleration of SEP may 
take place in the fiare region, solar corona and even in the interplanetary space. It should be kept 
i n mind that SEP them selves, while propagating, create conditions for the energy redistribution 



( Tvlka 



2001 



Led I2OO5I ). Therefore, the energy spectrum of SEP provides valuable information 



for study of solar and interplanetary plasma processes. Most energetic particles are usually (not 
always) accelerated during short time close to the energy explosion, thus the first particles arriving 
at the observer site keep more information about the primary acceleration process. Less energetic 
SEP may be accelerated during hours and even days. Conditions for particle acceleration in the 
corona and interplaneta ry space are determined by temporal and spatial evolution of a shock (e.g. 
Elhson fc R,amatvlll98d ;). 



Since the energy range of SEP extends over more than 5 orders of magnitude, and the SEP 
fluxes occupy more than 8 orders of magnitude, the SEP energy spectrum has never been measured 
by a single instrument. The bulk of the SEP observations has been performed on spacecraft and 
covers the particle energy below several hundreds MeV/n, whereas the relativistic solar particles 
have been observed by ground-based installations and remain less elucidated. This kind of obser- 
vations and associated modeling can also be used to issue radiation alerts to be used by manned 
and unmanned missions. A multi-detector study of the solar particle event occurring in 2005 and 
the correspon ding GLE69 was pe r forme d on-ground using NM and SREM radiation environment 
units in space iPapaioannou et al.l (1201 ll ) 



It is not yet clear if relativistic solar particles require a special scenario for their generation. 
Some authors believe that relativistic solar protons have two components - the fast and the delayed 
ones. T he fast component is initiall y accelerated in the flare region (e.g. IVashenyuk et al.l 120061 . 
2008al Jbl: iMcCracken Sz Moraall l2008l ) . probably via magnetic field reconnection, whereas another 
mechanism, p robably sho ck or stochastic acceleration, generates the second particle component. 
Bieber et al.l (l2004l . l2005l ) argue against two components and in favour of the CME-driven shock 
as a main accelerator of the whole SEP population. A wide variety in concomitant acceleration 
conditions supports the point that various in echanisms play the main role in different SEP events 



(Bombardieri et al. 



200. 



Bazilevskaval l2009l ) . Detailed studies of more solar particle spectra in 



wide energy range are needed to resolve this problem. 

Since ground-based instruments can only detect secondary cosmic rays, knowledge of re- 
sponse to the pr imary particles is ne c essary in order to f ind SEP fluxes in space. Some methods 



were developed (IShea fc Smartlll982l : I Cramp et al. 



1997 



Bieber et al. 2004: Plainaki et al. 



2007 



Vashenvuk et al.ll2006l ) to derive energy spectra and angular anisotropy of solar particles using the 
enhancements in the count rates of neutron monitors. Nevertheless, the particle fluxes derived from 
the neutron monitor data are model dependent. 

Before the PAMELA launch, direct measurements of the relativistic solar particles were not 
fulfilled. The aim of this paper is to present the results of the first direct measurements by a single 
instrument of the solar protons in the energy range from 80 MeV to several GeV and the helium 
nuclei from 75 MeV/n up to ~ 700 MeV/n during the SEP events of December 13 and 14, 2006. 
The absolute intensities and energy spectra are compared with the results of direct (GOES, ACE) 
and indirect (neutron monitors, IceTop) measurements. 



2. Instrument Description 

The PAMELA spectrometer (Figure [1]) consists of a number of higiily redundant detectors 
capable of particles detection trough the determination of their charge, mass, rigidity and velocity 
over a very wide energy range. Total weight of PAMELA is 470 kg; power consumption is 355 W. 



A more detailed descrip tion of the device and the data handling can be found in (jPicozza et al 
20071 : ICasolino et al.ll2006b| ]a). The core of the instrument is a spectrometer, constituted by per- 



manent mag net with a silicon microstrip tracker providing momentum and charge (with sign) 
information ( Adriani et alj|2003l ). The permanent magnet is 43.66 cm high and constructed of 5 
modules with high residual magnetic induction (0.43 T) provid ing an almost constant magnetic 



field value in the cavity. The tracking system (JRicciarini 



20071 ) is composed of 6 planes of high- 



precision Si microstrip detectors, positioned between the 5 magnetic modules of the tower, with 
uniform vertical spacing of 8.9 cm. Each silicon plane performs measurements of energy release 
and track position with a precision of about 3.0 fxm (X - bending view) and 11.5 /im (Y view). 



A scintillator system (JBarbarino et al.ll2003l ) provides trigger, charge and time of flight (TOF) 
information. It is composed of three double layers SI, S2, and S3, divided in various bars for a 
total of 48 scintillator paddles. It is also used to reject albedo particles which cross the detector 
from bottom to top. 

An anticoincidence system (CARD, CAT, and CAS) is used to reject spurious events in the 
off-line phase. 



A silicon-tungsten calorimeter (JBoezio et al.ll2002l ) is used to perform hadron/lepton separa- 
tion. It is composed of 44 silicon layers interleaved by tungsten planes for a total of 16.3 Xq 
(radiation lengths) and 0.6 Aj (nuclear interaction lengths). 

A shower tail catcher S4 and a neutron detector at the bottom of the apparatus are used to 
increase lepton/hadron separation. 

Around the detectors are housed the readout electronics, the interfaces with the CPU and all 
primary and secondary power supplies. All systems are redundant with the exception of the CPU 
which is more tolerant to failures. The system is enclosed in a pressurized container located on one 
side of the Resurs-DKl satellite. 



3. Data Selection and Analysis Criteria 

3.1. Geometrical Factor and Top of the Payload correction 

The Geometrical factor (Gf) of the detector has been evaluated defining a fiducial area con- 
sisting of a frame of 0.15 cm from the walls of the 13.1 x 16.1 cm? wide magnetic cavity. Only 
particles inside this fiducial area have been selected. This reduced volume ensures that all particles 



entering the magnetic cavity cross the scintillators and do not hit the magnet walls. The value of 
Gf = 19.93 cm?sr, constant within 1%, has been estimated with a numerical calculation and cross- 
checked with a Monte Carlo simulation where all the physical interaction processes are inactive, 
with the exclusion of particle bending in the magnetic field. The two values of Gf have been found 
in agreement within 0.1%. 

Interactions losses, due to local interactions where taken into account as a scale factor added to 
the Geometrical Factor. Protons and helium nuclei may be lost due to scattering and/or hadronic 
interactions in the 2 mm Al thick pressurized container in which PAMELA is housed or in the top 
of the detector. Correction factors amounts to ~ 6% for protons and ~ 12% for He due to the 
different cross section of the two species. 



3.2. Trigger system 

In low radiation regions, close to the geomagnetic equator and outside the South Atlantic 
Anomaly, where PAMELA crosses the trapped protons of the inner Van Allen belt, particles must 
cross at least one of the two layers of the three scintillator systems (51, 52, 53) to provide a valid 
trigger. In low cutoff regions or in the Van Allen belts, where particle rate is higher, 51 signal is 
not required to avoid random triggers due to the high number of particles. Given the low energy 
of solar particles compared to the galactic nuclei, in this analysis we have used particles selected in 
high latitude regions and thus taken with only (52 and 53) configuration. 



3.3. Live Time 

The live time tu^e of the apparatus is evaluated using the scintillator and trigger system. The 
counters for the live and dead time {tuve-, tdead) are cross-checked with the on-board time of the 
CPU, measuring the acquisition time {tacq = tu^e + tdead)-, to remove possible systematic errors due 
the counting method. The error associated with clock resolution is negligible compared to other 
sources of systematic errors. 



3.4. Event Selection 

3.4--1- Time of flight system selection 

In this analysis we have selected events that do not produce secondary particles in the 51 and 
52 scintillators and in the tracker, requiring a single fitted track within the spectrometer fiducial 
acceptance and a maximum of one paddle hit in 51 and 52 matching the extrapolated from the 
tracker trajectory. 



Particles interacting in the satellite can produce showers with different secondaries hitting 
the scintillator pads. These showers may produce random coincidences in the scintillators: they 
are rejected by means of anticoincidence and TOF cuts. We require the absence of hits in the 
anticoincidences around the empty area between 51 and 52 and around the top magnetic cavity 
(CARD and CAT respectively). No constraints on the anticoincidences around the magnet (CAS) 
and 53 have been put, since they are more often hit by backscattered secondaries produced in 
the calorimeter. The probability for such particles to hit CARD and CAT has been estimated 
with experimental data and cross-checked with Monte Carlo simulations. This efficiency has been 
included in the proton flux estimation. No constraints on the anticoincidence have been imposed 
for helium. 

The timing information of the TOF scintillator paddles along the extrapolated trajectory is 
used to evaluate the /3 of the particle. Albedo particles crossing the detector from bottom to top 
are discarded by requiring a positive (3. 



3.4-2. Proton and Helium identification 

Since average energy loss of a charged particle through matter follows Bethe Bloch formula, 
dE/dx oc Z"^ / j3'^ (neglecting logarithmic terms), the measurement of the average energy released 
in the tracker planes for a given event at a given rigidity can be used to identify different parti- 
cles. Proton and helium candidates have been selected requiring energy loss in the tracker planes 
compatible with Z=l and Z=2 nuclei. 

Cuts in the energy loss (dE/dx) vs. rigidity remove positrons, pions and particles with Z >2 
as shown in Figure [2l The bands in the Figure due to protons and helium nuclei which have energy 
loss in the tracker Z^ = 4 times protons, are identified. Events with small energy losses below 1 
GV are due to positrons, relativistic also at low rigidities and the background of secondary pions. 
The contribution from these two particles are negligble respect to protons above 1 GV. The black 
lines show the rigidity dependent cuts used to select the proton and helium samples. From the 
same figure the deuterium contribution at low rigidities can also be identified, resulting in a band 
with energy releases higher than protons. In this work we did not discriminated between isotopes. 

Using the redundancy of energy loss measurements in tracker and TOF, residual contamination 
of protons in the helium sample has been found to be below 0.5% in the energy range used for this 
work. 



3.4-3. Geomagnetic Selection 

The high inclination (70°) orbit of the Resurs-DKl satellite allows particles of different origin 
and nature to be studied. To separate the primary (solar and galactic) component from the reen- 



trant albedo component (particles produced in cosmic ray interactions with the atmosphere below 
the cutoff and propagating along Earth's magnetic field l ine), we evaluated the local geomagnetic 
cutoff G in the Stormer approximation ( Shea et al.l 119871 ) . The value of G = 14.9/L^ - valid for 
vertically incident particles - is estimated calculating the Mcllwain L shell with IGRF magnetic 
field model along the orbit (JNOAAll2011bl ). Particles were selected requiring R > 1.3 G to remove 
any effect due to directionality in the detector and Earth's penumbral regions. 



3.4-4- Tracker selection 

Particle rigidity is obtained fitting it's track in the spectrometer. In this analysis we selected 
events with a single track fully contained inside the fiducial acceptance. 

For each particle, the tracking system provides up to 12 position measurements (6 in the 
bending view), which are interpolated with a trajectory evaluated by integrating the equations of 
motion in the magnetic field. At low rigidies (below ~ 1 GV) we have corrected for the energy loss 
in the detector. 

The Maximum Detectable Rigidity (MDR) for a given detector is defined as the rigidity for 
which the relative error on the rigidity AR/R = 100% and varies between 200 GV and 1.5 TV, 
depending on delta-ray production and event's topology. In this analysis, therefore, the particle 
rigidity is well below the MDR and thus is no MDR requirements have been applied. 



3.4-5. Selection efficiencies and residual contamination 

The tracker efficiency has been measured selecting a sample of events which leave straight 
tracks in the calorimeter and do not interact hadronically. These tracks were propagated back 
through PAMELA acceptance and tracker efficiency has been evaluated. 

The total selection efficiencies for both species has been obtained as a product of Tracker, TOF 
(anticoincidence) efficiencies. The resulting efficiencies for proton and helium nu c lei aft er all above 
mentioned cuts are shown in Figure El For a detail discussion see lAdriani et al.l (120111 ) . 



Although we put strong requirements on Tracker, TOF and anticoincidence systems, never- 
theless there is residual contamination of secondary particles produced on the top of the payload 
entering the PAMELA acceptance window and passing selection cuts. The maximum contribu- 
tion to the background for protons comes from secondary single-charge particles (positrons and 
pions). In case of helium, rejection power of our selection cuts is enough to make residual con- 
tamination negligible. In order to estimate this contamination we carried out a 27r Monte Carlo 
simulation of the protons and helium nuclei hitting t he PAMELA pressurized container. Two d iffer- 
ent hadr onic interactions pa c kages, based on Fluka (JBattistoni et al.l 120071 : iFerrari et al.ll2005l ) and 
Geant 4 (jAllison et al.ll2006l : lAgostinelli et al.ll2003l ). have been employed to simulate these inter- 



actions. This background decreases with increasing energy and amounts to less than few percents 
at 1 GV. Flux attenuation was estimated with an energy dependent simulation and it is constant 
above several GV. 



4. December 13 and 14, 2006 Solar Particle Events 

4.1. General description 

The most significant SEP fluxes detected by PAMELA were started on December 6, 13 and 14, 
2006 and originated from the active region NOAA 10930. December 6 event originated at the East 
limb resulting in a gradual proton event reaching the Earth on December 7 and lasting until the 
events of December 13 and 14. Due to a scheduled maintenance procedure, no data were collected 
by PAMELA tracker during December 6 event, requiring a different analysis approach that will 
be discussed in a forthcoming paper. On December 13 200 6, at 0214 UT an X 3.4/4B solar flare 



occurred in the same active region NOAA 10930 (S06W23) (|N0AA - STPll2006l ). The intensity of 



the event was quite unusual for a solar minimum condition. This event produ ced also a full-hal o 



2011 



Coronal Mass Ejection (CME) with the sky plane projected speed of 1774 km/s (JLaboratory 

The forward shock of the CME reached Earth at 1438 UT on December 14, causing a Forbush 

decrease of galactic cosmic ray intensities which lasted for several days. 

The flare XI. 5 (S06W46) at 2107 UT on December 14 gave start to a new growth of particle 
intensity as recorded by PAMELA and other satellites. The maximum energy of protons was below 
1 GeV, therefore no ground level enhancement (CLE) was recorded. The corresponding CME had 
a velocity of 1042 km/s. 

Figure H] shows the solar particles intensity time profiles with various energies as observed 
by the PAMELA spectrometer, GOES-11 and neutron monitors. Ground-based neutron monitors 
started recording of SEP at 0248 UT, while the near-Earth satellites responded later to SEP of lower 
energies. First solar protons arrived at the Earth's orbit with anisotropic pitch-angular distribution 
which resulted in different intensity-time profiles at neutron monitors situated at location with the 
same geomagnetic cutoff Re but different asymptotic directions of incoming particles. A highly 
collimated bun ch of relativistic solar p articles was observed by the muon hodoscope URAGAN a t 



0256-0304 UT (JTimashkov et al.ll2008l ). Anisotropy vanished by 0430 UT (IVashenvuk et al.ll2008al ). 



The whole energy range of SEP detectable by the PAMELA spectrometer can only be taken in 
the high latitude regions, where the geomagnetic cutoff is lower than the rigidity of solar particles. 
Measurements were performed during 5 passages of PAMELA through the polar regions on December 
13. Data were missed from 1000 UT of December 13 till 0914 UT of December 14 because of onboard 
system reset of the satellite. 

The results of PAMELA measurement during the December 13 and 14, 2006 SEP events are 
shown in Figure [5] and [6] respectively. The quiet time galactic spectra were measured in November 



2006. The first (southern) polar passage of PAMELA started at 0318 UT and showed strongly 
flattening of p and He spectra at energies ~200 MeV/n. The helium spectrum in the initial phase 
of the December 13 event is almost similar in shape to the proton one, but flatter at energies 
below 200 MeV/n. In the event of December 13 the initial spectra were rather hard and changed 
significantly at least till 1000 UT. Further evolution of the proton and helium spectra was rather 
smooth even in spite of a Forbush effect in the middle of December 14, until arrival at ~2300 UT of 
the SEP accelerated in the flare at 2107 UT on December 14. The intensity of ~100 MeV protons 
firstly increased till -^0400 UT on December 13 and stayed at almost the same level for at least 4 
hours, while higher energy particles were rather rapidly removed from the Earth's vicinity. After 
that, the particle intensity decreased, but the spectrum shape was almost constant. Flux intensity 
did not decrease to the solar quiet levels because of the onset of December 14 event. During this 
event spectral form slightly flattened in the lower energy part thus remaining quasi-exponential. 
The highest proton energy for this event was around 500 MeV. 

We have examined temporal behavior of proton-to-helium ratio in the both solar events. For 
this analysis the rig idity spectra were evaluated with rigidity intervals chosen in accordance with 



Adriani et al.l (J201ll ). Figure [7] presents the proton-to- helium ratio obtained in selected time inter- 
vals. The galactic background was not subtracted and is shown as a black curves. It is seen in the 
left panel of Figure[7]that at the beginning of the December 13 event the ratio was almost the same 
as for GCR (note that SEP intensity was highest at that time). The subsequent measurements 
showed the enhancement of He relative to protons up to the late December 14 when the high rigid- 
ity He virtually disappeared after subtraction of the galactic flux. Such a behavior may to a first 
approximation be interpreted as a result of faster removing of solar protons from the Earth's orbit. 
Protons have higher velocity than He of the same rigidity therefore the diffusion coefficient of He at 
the same rigidity is smaller and protons leave faster. Further analysis needs a particle propagation 
modeling. Right panel of Figure [7] clearly shows that no solar helium with rigidity above 1 GV was 
observed in the event of December 14. 

No analysis directed to search of ^H and "^He has been performed yet. No electrons or positrons 
of solar origin were found by the PAMELA. However, it should be noted that the energy threshold 
of PAMELA is at ~ 50 MeV, therefore rather high for solar electrons. 



4.2. Initial phase of the December 13, 2006 event 

The first arriving particles are objects of special interest because they provide more information 
about mechanism and conditions of acceleration. As it was mentioned, solar particles at the initial 
phase of the December 13 event demonstrated highly anisotropic angular distribution. In such 
conditions the asymptotic directions of particles arriving at PAMELA should be taken into account. 

The geomagnetic field acts not only as rigidity analyzer but also as an angular analyzer of 
charged particles coming from space. Because of declining in the geomagnetic field a particle with 



10 



certain rigidity above the geomagnetic cutoff Re can arrive at a given site only from a definite 
direction outside the magnetosphere (asymptotic direction). Therefore, at any point of its orbit 
PAMELA accepted particles with certain rigidity from only one asymptotic direction. 

At the moment of the first particle arrival PAMELA was at the latitude of 40N moving to the 
equator. Therefore the low energy particles could not reach the PAMELA detector. In addition, 
viewing directions of PAMELA were changed rapidly while the spacecraft moved along the orbit. 
This can be seen in the upper panel of Figure [8] where the asymptotic directions of 7 GV and 10 
GV protons arriving at PAMELA are shown for the initial phase of th e SEP event. The a symptotic 
viewing directions of PAMELA has been calculated with the program of lGvozdevskv ( 200EJ ) using the 
back-tracing of solar particle tr ajectories and the mag netospheric model of lTsvganenkd (J2002l ). It 
is in reasonable agreement with lPlainaki et al.l (l2009bl ) obtained with the Tsyganenko 1989 model. 
The dir ection of the SEP anisot ropy axis, i.e. direction from which the bulk of SEP comes is also 
plotted (jVashenvuk et al.l l2008al ) . According to Figure [HI PAMELA viewing directions at the early 
phase of the SEP event would allow the solar protons in the range of 7-10 GV to be recorded by 
PAMELA at 0248-0250 UT as these particles arrived within a ~30 degree cone around the SEP 
flux axis. However, above 3 GeV the statistic is poor and we have found no increase of proton flux 
there. Changes in the viewing directions of PAMELA resulted in rapid receding of its favorable 
direction of acceptance from th e solar proton an i sotrop y axis. By 0256 UT (when a collimated 
particle bunch was observed bv iTimashkov et al.l (J2008l )) the geomagnetic cutoff of PAMELA was 
12 GV and only protons with rigidity between 15 and 17 GV could come from directions close to 
the anisotropy axis. 

The lower panel of Figure [8] shows the asymptotic directions of particles with different energies 
arriving at PAMELA in the beginning of the event and at the most optimal time of observation 
when the energy spectrum of particles with energies above ~80 MeV/n was measured. For com- 
parison, the asympt otic directions for the Apatity neutron monitor and for the IceTop installation 
(JAbbasi et al.l |2008J) are shown in Figure [H 



4.3. Comparison w^ith other data 

In Figure [9] the differential spectra of PAMELA with the galactic background subtracted are 
compared with other experiments. Each panel corresponds to different time period. 

In the low-energy r ange PAMELA particle fluxes can be compared with other spacecraft mea- 
surements. The GOES (JNQAAll2011al ) fluxes are in moderate agreement with PAMELA ones, with 
exception for the most energetic GOES band (160-500 MeV) where GOES fluxes are several times 
higher. The same is true for GOES helium spectra. This could be due to contamination of lower 
energy particles in these GOES channels. On the other hand, the PAMEL A helium fluxes are in a 
good agreement with measurements, performed by ACE/SIS (|ACEll201ll ) in the adjacent energy 
interval. 



11 



At high energy, our data can be compare d with the IceTop shower array (jAbbasi et alj 120081 ) 
and the NM network (|Vashenvuk et al.ll2008al ). The results of the IceTop shower array are plotte d 
in Figure [9^, b specially taken at the time intervals of PAMELA observations (jKuwabaral l2009l ) . 
The IceTop installation has detected secondary particles, generated by the solar energetic particles, 
i.e. the calculated yield functions have been used to reconstruct the primary SEP fluxes. The 
IceTop results are depicted only in the energy range over which this detector was actually seeing a 
significant number of particles. The SEP fluxes obtained by PAMELA and those obtained by IceTop 
is in relatively good agreement. However, at 0318-0329 UT (Figure Eh) the power-law shape of the 
IceTop spectrum does not match the PAMELA spectrum which is more exponential-like. At 0406- 
0420 UT (Figure [OJd) the spectra of PAMELA and that of IceTop are virtually coincident in the 
overlapping energy range. 

The energy spectra from the NM network have been derived from the early phas e of the event 
(before 0300 UT on December 13) tiU ~ 0500 UT on December 13 (|Balabinll2009l ). As already 
mentioned, PAMELA missed the early phase of the event and took the first measurement at 0318- 
0329 UT, during mildly anisotropic particle arrival (Figure [9^). A shape of the PAMELA proton 
spectrum is close to that derived from the neutron monitor network at energies between 700 MeV 
and 2 GeV. At lower energies the PAMELA spectrum is harder; at E > 2 GeV the uncertainty in 
the PAMELA spectrum is high because of small statistics. The absolute intensities of protons from 
PAMELA are somewhat lower than those derived from neutron monitors even at energies above 
1 GeV. To some extent, discrepancy between the PAMELA and neutron monitor spectra may be 
caused by the flux anisotropy. 

The GLE parame ters for this event have also been calculated using the NM-BANGLE model 
The direct comparison of the particle fluxes is difficult since these data in 



Plainaki et al. 



Plainaki et al. 



(2009a 



2009al ) are given only for integral energy spectra (see their Figure 3). An estima- 



tion shows that the resu lts of the NM-BANGLE model are very close to the spectra obtained by 



Vashenvuk et al.l (J2008al ) and therefore differ from the PAMELA results in similar manner. Figure 



m shows that during the first polar passage of PAMELA (0318-0329 UT) the Apatity and Barents- 
burg neutron monitors demonstrated different rate enhancements, i.e. the SEP fluxes were still 
anisotropic. The neutron monitor spectrum refers to the direction of the anisotropy axis where 
the proton intensity is maximal. The PAMELA viewing directions for the protons below 3 GV (see 
lower panel of Figure El) were not very close to the anisotropy axis, therefore the PAMELA fluxes 
may be lower than derived from neutron monitors. By the next PAMELA polar passage (0406-0420 
UT) discrepancy between the PAMELA and neutron monitor spectra increased (Figure [9)3) . Since 
the flux anisotropy becames lower, one should expect better consistency between the PAMELA and 
neutron monitor data. On the contrary, observations show discrepancy which persisted at around 
0500 UT (Figure inh) . Particle fluxes changed slowly at 0800-1000 UT, so in order to increase statis- 
tical significance, the PAMELA results were averaged over 0850-0944 UT (Figure. Ol). By this time 
the effect on the neutron monitor network was already too small for determining the solar proton 
spectrum. However, the PAMELA results were confirmed by an independent observation of the LPI 
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balloon-borne detector ([Stozlikov et al.ll2007l ) which measured solar proton energy spectrum in the 
stratosphere at 0946-1046 UT. Energy of solar protons is determined in the LPI experiment from 
absorption of protons in the air, so this is also a dire ct measurement of solar protons although in 
the limited energy interval ([Bazilevskava et al.ll20ld ). The good agreement between proton data 
of the two instruments in the energy range 90 MeV - 400 MeV is an additional check that our 
efficiency and systematics are under control. 



5. The energy spectrum fitting 



The observed spectral shape may give some indications on the acceleration mechanism (ICasolino et al 



20081 ). The proton energy spectrum measured in the early phase of the event, which is most close 
to the time of primary acceleration (0318-0349 UT), was fitted by the functions representative of 
various acceleration processes: 






^p = ARK2 (R/caT) 



1/2 



(1) 

(2) 
(3) 
(4) 



where $p is the proton fiux intensity, E is kinetic energy, R is magnetic rigidity, c is particle 
velocity, and K2 is a modified Bessel function of order 2, with a T as free parameter (a repre- 
senting an acceleration rate and T the escape time from the acceleration region). An exponen- 
tial in kinetic energv ([ll) or ri giditv ([2]) function is typical for simple models of D C acceleration 



(IVashenvuk et 



Krvmskii 



1977 



al. 



2006 



2008a 



b|), a po wer law is indicative of shock acceleration (lAxfordl Il981 



tion (jMcGuire &: von Rosenvinge 



Ellison fc Ramaty ligSSl ). and Bessel function is resulted from stochastic accelera- 



19841 ). The best fit parameters of approximation are presented 



in Table [H 



Function 


A, cm~'^sr-^s-^GV^^{GeV-^) 


Parameters 


xVndf 


A ■ exp{-E/Eo) 


(68.2 ±0.9) 


Eq, MeV 


(262 ± 2) 


3.0 


A ■ exp{-R/Ro) 


(179 ± 3) 


i?o, MV 


(366 ± 2) 


7.6 


A . R-i-m~Ro)/Ro^ ^g = iGV 


(13.1 ±0.1) 


7 


(2.70 ± 0.02) 


4.1 


5 


(1.99 ±0.06) 


A-R-{K2{R/c-a-T)Y/'^ 


(106 ± 2) 


aT 


(1.71 ±0.01) 


6.4 



Table 1:: The best fit parameters of the solar proton spectra in the energy interval 80 MeV 
GeV in the early phase of the event (0318 - 0349 UT). 



It should be mentioned that even if the resulting fits reproduce rather nicely the spectral 
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shape, x^ values are high for ah of them (see Tabled]), showing that these analytical formulas do 
not correctly describe the spectrum in the whole energy range. The best fit is obtained for the 
exponential in kinetic energy function ([1]) , however it fails to describe the highest and lowest energy 
tails where the observed spectrum is harder (in the -R > 2.5GV and R < 1.7 GV range). It proves 
that no single mechanism can describe the SEP energy spect rum in the wide energy interval. The 



previo us successful attempts of the SEP spectra fitting (e.g. ICramp et al.l 119971 : IVashenyuk et al. 
2008al ) were possible only for narrower energy range. 



Discussion and conclusion 



The PAMELA spectrometer was the first instrument which directly measured the relativistic 
SEP in the near Earth space. It is important since all previous such measurements were fulfilled 
with the ground-based installation and the derived SEP fluxes depended on the instrument response 
functions. The spectra of solar protons in the energy range 80 MeV - 3 GeV and helium 75 MeV/n 
- 1 GeV/n were measured during the first polar PAMELA passage at 0318-0329 UT. There is 
a good agreement between the protoris fiuxes measure d by PAMELA and those obtained by the 
IceTop installation (lAbbasi et al.ll2008l : lKuwabarall2009l ). K eeping in mind accuracy of estimation 



of absolute SEP intensities from the neutron monit or data (jVashenyuk et al.l 120061 : iPlainaki et al 



20071 : iBombardieri et al.l 120081 : iPlainaki et al.ll2009al ). reasonable agreement can be stated between 
the PAMELA and neutron monitor SEP fluxes. However, the PAMELA spectra are always harder in 
the low-energy interval probably indicating that neutron monitor yield functions are underestimated 
below ~700 MeV. During the second polar PAMELA passage the difference between the SEP fluxes 
taken from PAMELA and neutron monitors has become larger while agreement between PAMELA 
and IceTop remained very good. The correctness of the PAMELA observations was also confirmed 
by agreement with the direct SEP measurement on a balloon in the atmosphere. 

Evolution of intensity and spectral shape of relativistic sola r protons is often used to deri ve 
some information about SEP generation (Bombardieri et al.ll2008l : IVashenvuk et al.ll2006l . l2008a| ]bl: 
McCracken &: Moraalll2008l : lMoraal et al.ll2008l ). The indications were found that the first arriving 
relativistic particles would be accelerated in the flare region and have an exponential spectrum 
whereas the latter particles would be accelerated by a CME-driven shock and have a power-law 
spectrum. In this case an appropriate dynamics of the energy spectrum could be observed by 
PAMELA. Being at the low latitudes PAMELA missed the earliest anisotropic phase of the event 
of December 13. The results of the flrst PAMELA observation around 0320 UT on December 13 
conflrmed existence of a hard quasi-exponential spectrum expected from the magnetic reconnection 
in the flare region. However, in spite of changes in the SEP intensity a quasi-exponential spectrum 
persisted, actually till the advent of the newly generated protons of the December 14 event. It 
should be noted that in the energy interval relevant to the neutron monitor observations (> 500 
MeV) the spectrum could be fltted by a power law beginning from ~1000 UT, however at lower 
energies the spectrum was flatter. We did not flnd any spectral form which would flt the observed 
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spectra satisfactorily in the whole energy range and could prove a certain dominating mechanism 
of SEP acceleration. A quasi-exponential spectrum and fast temporal evolution of particle fluxes 
during several hours were present only in the event of December 13 when the relativistic protons 
were generated. In the event of December 14, without relativistic particles, a form of the solar 
proton energy spectra changed little and was almost power-law throughout the event. In addition, 
no Helium with energy above ~100 MeV/n was observed in the December 14 event. This may be 
indicative of special conditions leading to acceleration of particles up to relativistic energy. 
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Fig. 1. — : Scheme of the PAMELA instrument. SI, S2, S3 - TOF scintillator planes (each plane 
consists of 2 layers), CARD, CAT, CAS - anticoincidence system. Spectrometer - tracker surrounded 
with a permanent magnet, Silicon- Tungsten calorimeter, S4 - Shower tail catcher, neutron detector. 
Particles enter the detector from the top, cross the scintillators and are bent in the magnetic 
spectrometer before interacting with the calorimeter. 
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Fig. 2. — : The energy loss (dE/dx) vs rigidity in the PAMELA tracker. The black lines show the 
rigidity dependent cuts used to select proton and Helium samples. 
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Fig. 4. — : Intensity of solar protons (top panel) and helium (bottom panel) of various energies 
as a function of time. Top panel: Circles are the PAMELA protons in the ranges of 80-160 MeV 
(yellow), 160-500 MeV (blue), 500-1500 MeV (black), and 1.5-3.0 GeV. (green). Lines are the 
GOES 11 protons for 80-160 MeV(red) and 160-500 MeV (blue) ranges. The enhancement in the 
neutron monitor count rates in percents respect to the background before the event is given for 
Apatity (Rc=0.65 GV, violet line) and for Barentsburg (Rc=0.1 GV, brown line). Bottom panel: 
Squares are the PAMELA helium in the ranges of 75-125 MeV/n (black), 125-350 MeV/n (green), 
and 350-700 MeV/n (cyan). Black line is 75-125 MeV/n He from GOES 11. 



21 









Nov 2005 






13/0318- 13/0349 UT 


> 


: 


hM 


13/0433- 13/0459 UT 


« 10' 

0) 

"e 

S 10 

<A 

C 

f ^ 

Q. 

10-^ 


r 


-3^^ 


13/0818- 13/0917 UT 

14/1650- 14/2235 UT 








: 


^^^B^j^^ 




10' 


r 


q 


kkm^ 




10-' 


r 


< <i 1 


; Ipiy 





10-^ 



10 
Kinetic Energy (GeV) 



> 

O 10 



1 r 



I 10-^ 

"10-' 

10-^ 

10-' 

10 = 



Nov 2006 

13/0318- 13/0349 UT 
13/0349-1 3/0433 UT 
13/0433-1 3/0459 UT 
13/0818- 13/0917 UT 
14/1650- 14/2235 UT 




■ 13/0318 
13/0349 

■ 13/0433 

■ 13/0818 

■ 14/1650 



- 13/0349 UT xlO" 
- 13/0433 UT xlO' 
- 13/0459 UT xlO^ 

- 13/0917 UT xlO' 

- 14/2235 UT x10° 




10"' 10'^ 



1 2 3 4 

Kinetic Energy (GeV) 



^10^ 

8 10' 

<A 



10^ i- 



E 
o 

=^10' 

E 
3 

a 10 

1 

10-^ 

10' 



10" . 
5x10" 





13/0318 - 13/0349 UT x 10 

13/0349 - 13/0433 UT x 10 


• 




13/0433 - 13/0459 UT x 10= 

13/0818 - 13/0917 UT X 10^ 

14/1650 - 14/2235 UT x 10° 


r 

- .... 1 


. . . 1 





10"^ 2x10' 



Kinetic Energy (GeV/n) 



1 2 3 4 

Kinetic Energy (GeV/n) 



Fig. 5. — : Proton (top panels) and helium (bottom panels) spectra as measured by PAMELA for 
the SEP event of December 13, 2006. Left: Absolute fluxes, Right: After subtraction of galactic 
background. Fluxes are scaled to improve separation of curves. In the time slices of Dec 14, 
1650-2235 UT it is possible to see the Forbush decrease of galactic particles after the arrival of the 
CME. 
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Fig. 6. — : Proton and helium spectra as measured by PAMELA for the SEP event on 14 December 
2006. Left: Absolute fluxes, Right: After subtraction of galactic background. Fluxes are scaled to 
improve separation of curves. 
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Fig. 7. — : The proton-to-helium ratio vs. rigidity at selected time intervals during the solar event 
of December 13 (left panel) and December 14 (right panel). The ratios are scaled to clarify the 
results. The GCR background is not subtracted and the ratio observed before the solar events is 
given as a black line 
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Fig. 8. — : Upper panel: The asymptotic directions of 7 GV and 10 GV protons arriving at PAMELA 
in the initial phase of the SEP event at selected times of December 13, 2006 (given near the curves). 
Lower panel: Asymptotic directions of protons arriving at PAMELA, the Apatity and Barentsburg 
neutron monitors, the IceTop installation. The particle rigidity is marked near the corresponding 
direction. Red squares denote the direction of the SEP anisotropy axis. 
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Fig. 9. — : Energy spectra of solar particles measured by the PAMELA spectrometer (with galactic 
background subtracted) and by other experiments. In red there are the PAMELA protons, in cyan 
the PAMELA He, horizontal bars are the GOES and ACE SIS data, (a-c): green line is the spectrum 
derived from the neutron monitor network, blue line is the IceTop spectrum; (d): in blue there are 
the LPI balloon data. 



